Chaos on the conveyor belt 
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The dynamics of a spring-block train placed on a moving conveyor belt is investigated both 
by simple experiments and computer simulations. The first block is connected by spring to an 
external static point, and due to the dragging effect of the belt the blocks undergo complex stick- 
slip dynamics. A qualitative agreement with the experimental results can only be achieved by 
taking into account the spatial inhomogeneity of the friction force on the belt's surface, modeled as 
noise. As a function of the velocity of the conveyor belt and the noise strength, the system exhibits 
complex, self-organized critical, sometimes chaotic dynamics and phase transition-like behavior. 
Noise induced chaos and intermittency is also observed. Simulations suggest that the maximum 
complexity of the dynamical states is achieved for a relatively small number of blocks, around five. 

PACS numbers: 05.45.Pq 



I. INTRODUCTION 

Spring-block type systems are successfully used for 
modeling various complex systems which exhibit self- 
organization. Usually, the static friction coefficients ex- 
ceed the dynamical ones, in systems composed of two or 
more blocks connected by linear springs the coexistence 
of these friction types can lead to avalanches, nonlinear 
dynamics or Self- Organized Criticality (SOC) [ll-Q. 

Spring-block type models have broad interdisciplinary 
applications (for a recent review see [5|), and prove to be 
very useful in describing many complex phenomena in 
different areas of science. For example, they have been 
applied successfully to explain elements of the Portevin- 
Le Chatelier effect p , the fragmentation obtained in dry- 
ing granular materials in contact with a frictional sur- 
face [3-[2| , to understand the formation of self-organized 
nanostructures produced by capillary effects [10|, [ll|, to 
model magnetization processes and Barkhausen noise 
[l2i] , to describe glass fragmentation pJi], and even for 
modeling highway traffic [ij, Il5| . 

The first members of this model- family were presented 
in 1967 by R. Burridge and L. Knopoff ^] to explain 
the empirical Guttenberg-Richter law of the size distri- 
bution of earthquakes [13 • One of their original models 
(called BK model here) is composed of a chain of many 
blocks connected with springs and two planes, to model 
the sHding of tectonic plates (for a recent review see for 
example [18]). The blocks can slide with friction on the 
bottom surface and they are all connected by springs to 
the upper plane which is dragged with a constant veloc- 
ity. This model presents stick-slip dynamics, and the size 
distribution of the slipping events exhibits a power-law 
scaHng. 

Since the pioneering work of Carlson et al. y, |l9|, 
the dynamics of the BK model has been investigated 



< k k k k k 



:b 



•^ |zneda@phys.ubbcluj .rol 



FIG. 1. The scheme of the investigated system. Blocks of 
mass m connected by springs of spring constant k are placed 
on a conveyor belt that is moving with a constant velocity u. 



with several type of friction forces. In these studies both 
chaotic dynamics |2QM24J and phase transition-like phe- 
nomena [25| have been reported. 

The system used here has also been proposed by Bur- 
ridge and Knopoff [l6| , and later it was referred as " train 
model" [2^. As sketched in Fig. [H a spring-block chain 
is placed on a platform (conveyor belt) that moves with 
constant velocity. The first block is connected by a spring 
to a static external point. As a result, due to the dragging 
effect of the moving platform, the blocks will undergo 
complex stick-slip dynamics. 

The stick-slip dynamics of one block has been studied 
by analytical [27f| 3Q.] . numerical |28M3l| and experimen- 
tal methods [^, [32| . The results for a chain composed of 
several blocks are contradictory from the point of view 
of self-organized criticality [iH S HI, Hj. Undoubt- 
edly, chaotic dynamics has been found in such systems 
by many authors [22|, l35l - [37j . In these systems nonlinear- 
ity is introduced via friction forces, and several models 
have been studied from velocity- weakening [26l, 1271 , l35l , l36j 
to state-dependent [31] friction forces. It has been shown 
that with velocity- weakening friction forces (and a con- 
stant static friction), the case of two blocks is the simplest 
autonomous spring-block system exhibiting chaos [35|. 
It was also argued, that for systems composed of many 
blocks, chaotic dynamics and SOC can coexist |33j . 

It needs to be mentioned that a somewhat modified 
version of the train model has also been investigated, 
where in contrast with the former model, both the first 



and the last blocks are connected to static points. Ad- 
ditional springs are also introduced. This system also 
exhibits periodic and chaotic behavior [38|, [39| . 

Although the proposed problem has been thoroughly 
investigated for one, two and many blocks, we believe 
that a careful experimental and detailed simulation study 
for an intermediate number of blocks can reveal new 
and interesting dynamical complexity. Comparisons be- 
tween analytical, numerical and experimental results ex- 
ist mainly for the one block system. An exception is the 
work of Burridge and Knopoff p^] , where a chain of eight 
blocks was thoroughly investigated. The main purpose of 
their work was, however, to investigate the distribution 
of the potential energy released during avalanches, and 
they did not considered the problem from a dynamical 
systems view. Also, a general feature of previous numeri- 
cal studies is that the model parameters and assumptions 
are arbitrarily chosen, and usually the only argument for 
this is to make the calculations or the numerical code 
simpler. 

Here we consider an experimental setup composed of 
5 blocks, and computer simulations of a simple model 
for systems up to 10 blocks. In contrast to previous 
works, the model parameters are realistically estimated 
from experimental measurements. We also incorporate 
the surface inhomogeneities of the conveyor belt by us- 
ing a Coulomb type friction force varying stochastically 
about a mean. This feature appears to be essential since 
the temporal intermittency observed in the experiments 
cannot be recovered with only deterministic friction. A 
spring force with exponential cutoffs is used, to make the 
force profile more realistic and to avoid the collisions of 
the blocks with each other. 

To the best of our knowledge this intermediate system 
size has not yet been investigated experimentally. The 
choice of intermediate system sizes is motivated by the 
expectation that this is the range where the theory of 
collective phenomena and of dynamical systems overlap 
and thus tools borrowed both from statistical mechanics 
and chaotic dynamics become both apphcable. Despite 
the relatively small number of elements the system con- 
sists of, we surprisingly find a sharp phase transition-like 
behavior as a function of the conveyor belt velocity. In- 
terestingly, as the size of the system is increased this 
transition becomes smoother, a phenomenon that can be 
understood through the intermittency that was revealed 
in this region. Tuning the level of disorder to a cer- 
tain value, disorder induced phase transition-like behav- 
ior is also observed. We show that by using the Coulomb 
type friction forces this simple system exhibits a com- 
plex, chaotic dynamics accompanied by power-law type 
avalanche size distribution indicating SOC, noise induced 
intermittency and also noise induced chaos. 

The work is structured as follows. First the used exper- 
imental setup is presented and the results are described 
for a chain of five blocks (Sec. |II]). In Sec. Illll the model is 
detailed. Simulation results (Sec. lIVp without and with 
a disorder in the friction force are discussed. Finally, the 




FIG. 2. (color online) Snapshot of the experimental setup. 
The spring-block system is placed on the conveyor belt of a 
treadmill. 



influence of the number of blocks (up to TV = 10) on 
the observed dynamics is investigated (Sec. |Vl) and final 
conclusions are drawn (Sec. IVip . 



II. EXPERIMENTS 

The experimental setup is sketched in Fig. [H The 
chain is built by TV = 5 black wooden blocks of mass m = 
115.8 g, and of dimensions 4 cmx8 cmx4 cm. The blocks 
are connected by steel springs of rest length / = 7 cm and 
spring constant of k = 19.8 N/m. The chain is placed on 
the conveyor belt of a treadmill (running machine) with 
adjustable speed. A digital camera is placed above the 
chain to record the dynamics of the blocks with 24 fps 
(for a snapshot see Fig. [2j). In order to make the last 
block clearly distinguishable from the others, and from 
the black colored conveyor belt, the top of this block is 
colored white. 

The video recordings are converted by a threshold op- 
eration into black and white image sequences as shown 
in Fig. m Then, the position of the last block is detected 
on these images. The length unit on the image can be 
determined using the image of a tape measure, which is 
placed next to the chain. Accordingly, the length xs{t) 
of the chain as a function of time can be obtained auto- 
matically with a simple processing program. 

The experiments were carried out with two values of 
the belt's speed: u = 0.22 m/s and u = 0.28 m/s. 
The measured average values of friction forces in the 
experiment are Fgto = 1.98 N for static friction and 
F/eo = 0.89 N for kinetic friction. Thus, the ratio of 
the two types of friction forces is fs = Fk^/FstQ = 0.45. 

In both experiments the system is initialized with 
blocks sticking to the belt and with undeformed springs. 
Initially, the conveyor belt is at rest. After it is set in 
motion, length oscillations of very large amplitude occur 
on a time- interval up to 100 s. In order to ensure that 
these initial transients are not considered and a kind of 
steady state has set in, a time-interval of ttrans = 200 s 
is discarded in the data processing. The recorded data 
for the position of the last block x^{t) as a function of 
time is then analyzed. For the velocity u = 0.28 m/s, 
the recorded data reveal (Fig. [3]) two qualitatively differ- 




FIG. 3. A processed black and white image sequence of the 
center of mass coordinate xs of the last (5th) block. 



ent temporal behaviors. In the domain characterized by 
small amplitudes in X5(t), the dynamics of the system is 
nearly periodic, while in the domain of large amplitudes 
the last block exhibits a chaotic- looking behavior. We 
recall again, that the initial periodic-like behavior in this 
figure is not a result of the lasting transients, since be- 
fore this, several chaotic-like regions were already present 
during the large amplitude bursts following the start of 
the conveyor belt. 

The differences in the dynamics for these domains can 
be better understood from the Fourier Transforms (FT). 
We have thus computed separately the FT for the two 
well distinguishable regions (see Fig. |4]). It is clearly ob- 
servable that for the small amplitude interval the power 
spectrum S of the FT has peaks at equal distances sug- 
gesting periodic dynamics, while in the chaotic region the 
power spectrum presents a quasi-continuous distribution. 

As can be seen in Fig. [5l for ii = 0.22 m/s the dynam- 
ics of the system is intermittent. In the first 125 seconds 
there are large "avalanches" that result in a simultane- 
ous movement of blocks (they slip and stick to the belt 
together). In this regime, the length of the chain evolves 
in time with fluctuations of large amplitudes. Then, 
without any external influence, the dynamical behavior 
changes abruptly, and for more than one minute the sys- 
tem remains in a nearly steady state, when all the blocks 
are slipping continuously relative to the belt. After this 
time interval, the behavior of the system looks chaotic 
again, of the type it was at the beginning of the recorded 
dynamics. The power-spectrum of the FT for the whole 
plotted interval is similar to the one observed for the 
chaotic regime for u = 0.28 m/s. 



III. A SIMPLE MODEL 

The model contains the same elements (blocks and 
springs) as the experimental setup. The main challenge 
in the modeling effort is however the quantification of 
the friction and spring forces and the numerical integra- 
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FIG. 4. (color online) Experimental result for the chain length 
X5(t) as a function of time when the velocity of the conveyor 
belt is 1^ = 0.28 m/s. The power spectrum S of the FT is 
computed separately for the two clearly distinguishable re- 
gio 
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FIG. 5. (color online) Experimental result for the chain length 
X5(t) with u — 0.22 m/s. The power spectrum S of the FT 
for x^{t) is presented in the inset. 



tion of the equations of motion. Let us consider a chain 
formed by N identical blocks of mass m. The blocks are 
connected by identical springs of linear spring constant k 
and undeformed length /. As Fig. [1] indicates the motion 
of the blocks of the model is one dimensional and the 
transverse oscillations present in the experiment are thus 
neglected. 

In the model, dimensionless units are taken. We have 
chosen these units so that m = 1, /c = 1, and / = 50. 
The value for / was chosen for the sake of an easier 
graphical visualization (spring length corresponding to 
50 pixels). Being motivated by the experiment, we can- 
not avoid using a nonzero rest length for the springs. In 
order to have these dimensionless values to correspond 
to the experimental situation, the units are chosen as: 
[m] = 0.1158 kg for mass, [k] = 19.8 N/m for spring 
constant, and [/] = 1.4 • 10~^ m for length. The units 
of the other quantities follow from dimensional consid- 
erations . The time, velocity and force units are thus, 
[t] = y^[m]/[k] = 0.0765 s, [u] = [l]/[t] = 0.0183 m/s. 







■ \ 


ictcai spring 




■ 


"■^j 



150 

100 

50 

Fe(Al) 

-50 

-100 

-150 

-40 -20 20 40 60 80 100 

Al 

FIG. 6. The used nonlinear spring force profile, in comparison 
with the ideal spring force profile. 



and [F] = [k] • [/] = 0.0277 N, respectively. 

The equation of motion for the ith block of the chain 
is written as 

X, = Fe{Al_) - Fe(A/+) + Ff[Vr,,Fe{Al_) - i^e(A/+)], 

(1) 

where A/_ = Xi — Xi-i — 50 and A/+ = x^+i — Xi — 50, 
respectively, and v^ is the relative velocity of the block to 
the conveyor belt. The elastic force, Fg, and the friction 
force, Fj, are defined below. 

The elastic force F^ of any spring is linear, up to a cer- 
tain deformation value, Al^ax- For higher deformations, 
this dependence is assumed to become exponential, with 
an exponent bigger (in modulus) for negative deforma- 
tions (see Fig. [6]). Accordingly, 
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(2) 
where we have chosen Al^ax = 20, bi = 0.2 for A/ < 
and 62 = 0-01 for Al ^ 0. With these choices the model 
is more realistic since the nonlinearity of spring forces 
is taken into account and collisions between blocks are 
avoided. At the same time the choice of parameters bi 
and 62 is somewhat ad-hoc and we have not carried out a 
systematic change of them, other than the fact that the 
average chain lengths are thus comparable with those 
seen in the experiment with an error less than 10%. 

According to a review of experimental results pre- 
sented in Ref. [32|, "the classical friction law, where the 
friction force is proportional to the load, will only ex- 
ist in an average sense" and it is important to take into 
account that in many cases the friction is determined 
by surface asperities, which means that "deterministic 



friction- velocity relations at best only exist in an average 
sense" . In |32| , the authors also argue that in the case of 
stick-slip dynamics on a plane surface, the friction forces 
may have normal distributions. 

In the present paper Coulomb's law of friction [40| is 
used with a noisy extension reflecting the surface irreg- 
ularities. Both the static and kinetic friction forces are 
independent of velocity modulus. A block remains in a 
stick state until the resultant external force F^x exceeds 
the value of the static friction force, Fst- For higher ex- 
ternal force values the block starts to slide in the presence 
of the kinetic friction force F/.. We assume that, the ra- 
tio of the static and kinetic friction forces Fk/Fgt = fs is 
constant. The friction force Ff acting on the block de- 
pends both on sgn (vr), where sgn is the signum function 
and Vr is the block's speed relative to the conveyor belt, 
and on the value of the resultant external forces F^x act- 
ing on it. In our ID setup, the friction force orientation 
is given only by its sign: 
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(3) 
where Vr = v — u and v is the velocity of the block relative 
to the laboratory frame. 

Since the surface of the conveyor belt is not perfectly 
smooth, the already mentioned argument of Ref. [32[ 
is implemented by using randomly distributed friction 
forces. This is incorporated into the model by randomly 
generating a new static friction force value for every dif- 
ferent position of the block relative to the belt. These 
random friction force values are generated according to a 
normal distribution with a fixed mean Fgt^ and standard 
deviation: a. This standard deviation (together with the 
integration time-step, which is however the same for all 
simulations), will quantify the amount of disorder intro- 
duced in the model. Although the origin of the disorder 
is the spatial inhomogeneity of the belt's surface, from 
the point of view of the dynamics the fluctuations of the 
friction force appear as temporal noise. We shall thus 
call (7 as the noise strength. 

The kinetic friction force value is always automatically 
updated together with the static friction force value, us- 
ing their fixed ratio fg. As a result, both the kinetic 
and static friction forces will fluctuate in time during the 
sliding of the particular block. In Fig. [7] we illustrate 
the variation of the friction force as a function of the re- 
sultant external force F^x for one stick-slip event. Such 
type of friction forces also proved to be useful for model- 
ing highway traffic [5|, [ij] . 

In order to use the same friction force value as in the 
experiments, in the dimensionless units the average value 
of the static friction force is {Fgt) = FgtQ = 71.4. The 
ratio of the two types of friction forces is also chosen in 
agreement with the experiments as fs = 0.45. 

Computer simulations of the model start from initial 
conditions similar to the experimental ones. The blocks 
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FIG. 7. The stochastic version of the Coulomb friction model 
for one stick-slip event. Thick dashed line represents the situ- 
ation when the block is stuck to the belt (the resultant exter- 
nal force Fex is smaller than the maximal static friction force 
value Fst)- Continuous line corresponds to the situation when 
the block is sliding and indicates the fluctuation of the kinetic 
friction force Fk = fsFst during such an event. Dashed gray 
arrows indicate the sudden change of the friction force, when 
the block starts to slip. 



are placed on the belt with undeformed springs between 
them. At this initial setup the blocks are stuck to the belt 
which means that their initial velocity in the laboratory 
reference frame is equal to the conveyor belt's constant 
velocity u. Computer graphics is also used to visualize 
in real-time the dynamics of the blocks. 

The numerical methods used to solve Newton equa- 
tions ([!]) are presented in the Appendix. The time-step 
of integration is fixed to dt = 0.001. 

In order to characterize statistically the dynamics of 
the chain, a parameter measuring the fluctuation of the 
chain length in dynamical equilibrium is introduced. The 
chain's length is determined by the position of the last 
block in the row. The disorder parameter r is defined as 
the standard deviation of the coordinate xn compared 
to the time averaged length (xn) of the chain: 



^J{X%) - {xnY 

{xn) 



(4) 



This quantity takes on large values for large fluc- 
tuations which explains the term disorder parameter. 
Among other relevant dynamical measures, this disorder 
parameter is investigated as a function of the belt veloc- 
ity u (Section IIV A[ ). the noise level a (Section IIVBI ). 
and the number of blocks N (Section IVl). 



IV. SIMULATION RESULTS 
A. Results without noise 

In the deterministic case, as the conveyor belt is 
started, the whole system moves together with the belt 
until the first block starts to slip. This slipping moment 
is determined simply by the value of the static friction 
force Fst and the belt velocity u. The block sticks again 
to the belt when its relative velocity Vr becomes zero. 
After that the process starts again. For small belt ve- 
locities this kind of behavior characterizes all the blocks 
in a self-organized manner. Specifically, the blocks are 
slipping together and produce "avalanches" of different 
sizes. Therefore, the length of the chain, defined by the 
position of the last block xn fluctuates as a function of 
time with largely varying amplitudes, as shown in Fig. 
Ela. 

Furthermore, as shown in the inset of Fig. [8la. the 
Fourier transform of the time series of xn exhibits a 
nearly power-law behavior with an exponent —1.21, sug- 
gesting a 1// type stochastic behavior in the long time 
dynamics and a critically self-organized state. The rela- 
tive velocity of the last block Vr^ is zero, when the block 
is stuck to the belt. As the chain starts to slip, this 
relative velocity is rapidly increasing in absolute value 
as indicated in Fig. [8lb and it shows large fluctuations 
according to the inset of Fig. [8lb. 

For a quantitative measure of avalanches, the slip size 
distribution for the last block TV = 5 is computed and 
shown in Fig. [9la. The slip size Axg is defined as the 
difference of the coordinates where the block starts to 
slip and when it stops relative to the belt. The obtained 
power law behavior confirms the presence of SOC, which 
appears here together with a chaotic dynamics. The dis- 
tribution of the energies dissipated during avalanches is 
plotted in Fig. [9lb. The results indicates again a scaling 
(a sign for SOC-like behavior), although the exponent 
a = —0.87 is different from the one obtained in the sem- 
inal work of Burridge and Knopoff [16]. This difference 
is a natural consequence of the different nature and size 
of the problems. 

The existence of chaos is illustrated by Fig. [TOl where 
the natural distribution on a Poincare section is pre- 
sented. This Poincare section is projected onto the phase 
plane (xat, Vr^^) of the last block, and obtained as the in- 
tersections of the system's trajectory in the phase space 
with the plane defined by X2 = 304.3, considering only 
uni-directional crossings from right to left (±2 < 0). For 
a better view, the high probability states along the line 
^VN — (where the block is stuck to the belt) are not 
shown. In the full system of 10 degrees of freedom the 
Poincare section is 9-dimensional. Fig. [TOl shows a pro- 
jection of this on a plane. It is surprising that the natural 
distribution on this plane is similar in appearance (both 
in the shape of the support and of the rather irregular 
distribution) to that of the chaotic attractor of a driven 
one-dimensional system. This indicates that the chaotic 
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FIG. 8. (color online) Computer simulation results for u — 
1, (T — 0, and N = 5. a. The position of the last block 
xn as a function of time. Inset: the power spectrum of the 
corresponding FT. The fitted line corresponds to a power law: 
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b. The relative velocity Vrp^ of the last 



block as a function of time, 
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FIG. 9. a. Probability density function for the size, Axs, of 
the avalanches. A power- law fit with an exponent a = —1.47, 
suggests a SOC-like behavior, b. Probability density function 
for the energy A^ dissipated during the avalanches. A power- 
law fit with exponent a = —0.87 is shown. The simulation 
parameters are u = 1, a = 0, and A^ = 5. 



attractor of the full chain is rather low-dimensional. 

For higher values of the belt velocity, there are other 
possible scenarios as well. For instance, foTU = 7 the 
system exhibits an asymptotically periodic stick-slip dy- 
namics, but for u = 15 the behavior is aperiodic again. 
For the velocity interval IS < u < 50 permanent chaotic 
dynamics never occurs for the investigated parameters 
(see e.g. Fig. [11] for u = 20) and the last block exhibits a 
continuous slip dynamics after a certain transient time. 
It can be seen in this figure that the system undergoes 
a transient chaotic behavior with a stick-slip dynamics, 
before reaching a periodic attractor. 

The bifurcation diagram plotted in the top panel of 
Fig. [T2l summarizes all the cases described above. Here, 
the velocity of the last block vn is plotted as a function 
of the driving velocity u in instances when the block goes 
through a fixed {xn =615) position from the right. For 
every value of the control parameter u the simulation is 
restarted from the declared initial position. To construct 
this plot, the computer simulations were run up to t = 
2 X 10^ time units, discarding a transient time of ttrans = 
10^. 

Interestingly, the disorder parameter dH allows for a 
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FIG. 10. (color online) Natural distribution on the Poincare 
section in the phase plane of the last block (i^ = 1, a = 0, 
A/" = 5). The Poincare section is taken when the X2 = 304.3 
plane is crossed from right to left. The plot is made only for 
Vrj^ < (the stuck states with Vr^ = are not considered for 
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FIG. 11. a. The position of the last block as a function of time 
for u = 20. The inset presents a zoom on the indicated time- 
interval, b. The relative velocity Vrj^ of the last block. The 
last block exhibits a chaotic transient dynamics up to about 
t — 4000, followed by a periodic behavior with continuous slip. 
The inset presents a zoom on the indicated time-interval. The 
other parameters are N — 5 and a — 0. 
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FIG. 12. Top panel presents the bifurcation diagram for the 
velocity vn of the last block captured when xn = 615 and 
xn < 0, as a function of u. Bottom panel illustrates the 
disorder parameter r as a function of i^ (a = and A^ = 5). 
The inset presents the average lifetime r of transients for the 
18 < 1^ < 35 interval. The missing r values in the interval 
1^ < 18 indicate that the lifetime of the transients is infinite, 
chaos is permanent. 



namics is permanent, while for i^ > 18 it has a transient 
character. These transients were neglected in comput- 
ing the disorder parameter, since the discarded transient 
time {Urans = 10^) is larger by more than two orders of 
magnitude than the average lifetime r of the transients. 



difference to be made between these different dynamical 
behaviors. Results for r{u) are plotted on the bottom 
panel of Fig. [121 As expected, for asymptotic chaos 
the disorder parameter has a high value, and in the case 
of the periodic or quasi-periodic dynamics, it is much 
smaller. In the transition region (15 < i^ < 18) the value 
of r falls from about 0.13 to 0.02. This jump of nearly one 
order of magnitude indicates a relatively sharp dynamical 
phase transition-like behavior. The critical speed Uc^ de- 
fined as the midpoint of the transition region, which the 
phase transition-like behavior can be associated with, is 
Uc = 16.5. (Note that in the interval 20 < ix < 22 there 
are three outlier points in the lower panel of Fig. [T2j 
These belong however to periodic attractors of large am- 
plitudes rather than to chaotic cases.) 

In the inset of Figure [121 the average lifetime r of the 
chaotic transients preceding the periodic behavior is pre- 
sented. This is measured using the method described 
in Ref [41|. Looking at the simulation data, it can be 
seen that in the phase transition region the lifetime of 
transients grows considerable, a kind of critical slowing 
down can be observed. It can thus be concluded, that for 
u ^15 (except for the periodic windows) the chaotic dy- 



B. Results in the presence of noise 

As described earlier, the two friction forces are linked 
by the relation F^ = fsFst^ and the values of Fgt are 
randomly updated for each new position of the blocks on 
the conveyor belt. The distribution of Fgt has a standard 
deviation a and a mean value Fsto - 

First, a relatively low level of noise <j = 1 is considered 
(see Fig. [7|). In this case, the phase transition-like behav- 
ior remains almost unchanged. This is clearly visible in 
the r{u) plots of Fig. [131 The periodic windows present 
for (7 = (top panel of Fig. [12]) disappear, however, and 
consequently the disorder parameter has less fluctuations 
in the chaotic regime. The noisy bifurcation diagrams of 
Fig. [TH also illustrate the disappearance of the periodic 
windows as a increases. 

For belt velocity u = 7 and small noise levels, the 
asymptotic periodic dynamics is reached after a tran- 
sient chaotic behavior. This is possible if a non- attracting 
chaotic set (chaotic saddle) and periodic attractors coex- 
ist (see for example [42|). Beyond a critical noise level 
(Jc, the chaotic transients turn into a permanent chaotic 
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FIG. 13. (color online) The disorder parameter r as a function 
of the conveyor belt's velocity u for different noise levels a. 
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FIG. 14. (color online) The bifurcation diagram, as defined 
in Fig. I12[ for different noise levels. The inset presents the 
disorder parameter as a function of a for i^ = 7 (a periodic 
window for a = 0). 



dynamics. This is called noise induced chaos |4l|, l43l - 
|46J . Based on the results plotted in the inset of Fig. [M] 
the critical noise strength Cc necessary to obtain noise 
induced chaos is estimated as <Jc = 0.75. 

By increasing the noise level a in the friction force, 
the critical speed Uc{o-)^ for which a phase transition-like 
behavior occurs, increases with a as suggested by Fig. 

M 

Turning the problem around, in Fig. [15] we plot the 
disorder parameter r as a function of the noise level a for 
intermediate velocities u ^ 20. In the whole range a noise 
induced phase transition-like behavior emerges. Moving 
away from u = 20, the sharp transition becomes more 
smooth and the transition point in a is shifted toward 
higher noise levels. A possible explanation of this effect 
is that for higher velocities a higher noise level is needed 
to kick the system out from the basin of attraction of the 
periodic attractor. 

The model is able to reproduce the intermittent be- 
havior observed in the experiment and presented in Fig. 
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FIG. 15. (color online) The disorder parameter r as a function 
of a for different values of u ^ 20. The a = 5 line is only to 
guide the eye. 



21 For example, considering a = 2.2, and belt velocities 
in the phase transition-like region, the system exhibits 
an intermittent dynamics (Fig. [T6|) . Since such inter- 
mittency is not present without noise, this phenomena is 
called noise induced intermittency |42|] in the chaos lit- 
erature. This can happen again only if a chaotic saddle 
exists in the deterministic system [42|. 

We note that the critical speed Uc{o- = 2.2) is found 
to be 19 which corresponds in dimensional units to Uc = 
0.35 m/s. Experimentally we have found intermittency 
dit u = 0.22 m/s and u = 0.28 m/s which should belong 
to velocities close to the critical value. In view of the sim- 
plicity of the model this order of magnitude agreement is 
satisfactory. 

The disorder parameter can be computed also for the 
experimental time series of x^{t) of Figs. |4]and[5l The 
results are: r = 0.072 and r = 0.052 for u = 0.22 m/s and 
u = 0.28 m/s, respectively. The values corresponding 
in the model to these dimensional velocities are: r = 
0.132 and r = 0.134, respectively, again in an order of 
magnitude agreement with the experiments. 

Another facet of the effect of noise is shown in Fig. 
[TTl where the natural distribution is projected to the bi- 
furcation diagram defined in Sec. IIV A[ The consecutive 
graphs are results for increasing noise levels. From these 
plots we learn that for <j = 1 all the periodic windows in 
the < 1^ ^ 15 interval disappear, but the phase tran- 
sition point hardly changes. As a is increased further, 
the sharp phase transition-like behavior transforms into 
a smoother one (see Fig. [TJlc). For a = 10 the dynamics 
of the system is dominated by noise (Fig. [TTld). 

The results regarding the phase transition-like behav- 
ior are summarized in Fig. [18] by a detailed map of 
the parameter space {u^a). Figures [13] and [15] corre- 
spond to sections of this map along the horizontal and 
vertical axes, respectively. The range where the sys- 
tem exhibits noise induced intermittency is enclosed by 
white dots. The fact that this range does not reach the 
a = line shows, that intermittency cannot be recovered 
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time exhibits a similar intermittent behavior as found in the 
experiments. The inset shows the power-spectrum of the FT. 
iV = 5, 1^ = 17.5, and a = 2.2. 



periodic time intervals appear at increased chain veloci- 
ties. Therefore, the sharp transition becomes smoother. 
For increasing system sizes, the previously observed noise 
induced phase transition-like behavior becomes also less 
evident, as illustrated by Fig. [211 fo^ driving velocities 
above the transition region. 

Bifurcation diagrams generated for different system 
sizes are shown in Fig. [221 These results augment 
those obtained from disorder parameter. They indicate 
that for TV < 5, the system has only periodic or quasi- 
periodic attractors. The behavior of the disorder param- 
eter does not suggest any phase transition-like behavior 
here. Dominant chaotic regimes appear for N ^ 5. As 
the size of the system is increased, the number of periodic 
windows decreases, and the transition region is enlarged. 
Also, for larger systems the effect of noise is less pro- 
nounced. 



with purely deterministic friction forces. As the level of 
noise is increased, this region becomes abruptly wider 
for cr > 4. This can be interpreted as the critical noise 
strength below which the phase transition-like region re- 
mains sharp, and if this is exceeded, the transition be- 
comes smoother. In the view of all these observations 
we conclude that the noise strength corresponding to the 
experiments is about a = 1 — 2. 



V. FINITE SIZE EFFECTS 

We can now briefly discuss the influence of the sys- 
tem size N on the observed phase transition-like phe- 
nomena. We know from statistical physics that the ob- 
served phase transitions become sharper as the system 
size is increased. Although the investigated system is a 
non-equilibrium dynamical system and not a system in 
thermal equilibrium, one might expect that the phase 
transition-like behavior sharpens for larger system sizes. 
In our system, however, the opposite happens. As shown 
in Fig. [19] for increasing system sizes, the sharp phase 
transition-like behavior is transformed into a smoother 
and smoother one, the transitional region broadens with 
N. 

This observation can be explained through a quanti- 
tative change in the systems dynamics in the presence 
of noise. As we have seen before, the dynamics changes 
from chaotic to periodic or quasi-periodic via an intermit- 
tent behavior. For A^ = 10, intermittency occurs even 
without noise. This is nicely illustrated by the XN{t) 
graphs for different u values in Fig. [20l The intermit- 
tency found in the transition interval 27 < u < 40 ex- 
plains why the sharp phase transition-like behavior disap- 
pears for larger systems. As the disorder parameter mea- 
sures the fluctuations of the chain length, in the intermit- 
tent region smaller disorder parameters will be observed. 
Fig. [20] indicates that larger and larger periodic or quasi- 



VI. DISCUSSION 

The dynamics of a simple spring-block chain placed on 
a running conveyor belt was investigated both by sim- 
ple experiments and through computer simulations. De- 
spite its simplicity, the dynamics of the system proved to 
be quite complex exhibiting- chaotic, periodic or quasi- 
periodic behavior as a function of the conveyor belt's 
velocity and the amount of noise in the friction forces. 
The experiments and computer simulations indicate that 
the transition from chaotic to periodic or quasi-periodic 
dynamics is typically realized through an intermittent 
chaotic state. 

In the chaotic regime, the avalanche-size distribution 
function shows a scale- free nature, indicating also the 
presence of SOC, along with a 1// type of stochasticity 
in the chain length. Another aspect of the systems com- 
plexity is the diversity of the possible dynamical states 
reflected by the convoluted graph of the bifurcation dia- 
gram, and a phase transition-like behavior in the disorder 
parameter. The presence of noise adds to this picture fas- 
cinating phenomena, like noise induced chaos, noise in- 
duced intermittency and noise induced phase transition. 

Interestingly, the computer simulations suggest that 
the maximal complexity of the dynamical states is 
achieved for a relatively small number of blocks in the 
vicinity of iV = 5 for <j < 2. For N < 5 there is practi- 
cally no chaos in the system. With A^ = 5 the observed 
transition is the sharpest one, and for A" > 5 it becomes 
smoother in a somewhat counter intuitive way. It is re- 
markable that this collective behavior finds its explana- 
tion in terms of dynamical system properties. 

Finally we would like to draw the attention to the fact 
that in several recent experimental studies [481, E9| per- 
formed on metallic alloys, it has been observed that just 
before the onset of the plastic instability the local strain 
rate as a function of the corresponding force response 
shows a multi periodic behavior. This is somehow sim- 
ilar with the graphs shown in Fig. [Til This observa- 
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FIG. 17. (color online) The natural distribution on the plane defined by u and vn for increasing noise levels (N=5). a. a = 0; 
b. cr = 1; c. cr = 5; d. a = 10. 



10 
8 
6 

4 

2 




LL.^ 



r 

i-i 0.16 

0.14 

0.12 

0.1 

0.08 

0.06 

0.04 

0.02 





10 15 20 25 30 35 40 45 50 
u 



FIG. 18. (color online) Disorder parameter r over the pa- 
rameter plane {u, a) for A^ = 5 blocks illustrated with color- 
coding. White dots indicate the boundaries of the parameter 
range where noise induced intermittency can occur. 



tion brings us back to the possibility of modeling the 
Portevin- Le Chatelier effect with a spring-block system 
[6|. Indeed, if the plastic deformation is macroscopically 
uniform in the absence of plastic instability, one can as- 



0.2 



0.15 



r 0.1 



0.05 



' 1 • 

o N=5 
A N=7 
□ N=10 




FIG. 19. (color online) The disorder parameter r as a function 
of the belt velocity u for different systems sizes N with the 
a = 1. The "transition regions" corresponding to different 
system sizes are: 15 < u < 18 for A/" = 5, 23 < n < 32 for 
AT = 7, and 27 < IX < 40 for N = 10. 



sume that it is governed by only a small number of col- 
lective degrees of freedom, similarly with the dynamics of 
the spring-block chain considered here. In such view the 
studied model and the obtained results can gain interest 
in the field of material science as well. 
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Appendix: Numerical method 

We briefly describe here the method used to integrate 
the Newton-equations ([1]) and to handle the discontinu- 
ous stick-slip dynamics of the blocks. If a block is stuck 
to the conveyor belt, then it will move together with it 
at constant velocity u. Therefore, the position of the ith 
block relative to the ground is calculated with the simple 



Xi{t + dt) = Xi{t) -\- u ■ dt^ 



(A.1) 



equation. When a block is slipping relative to the belt, 
the basic Verlet method 

Xiit^dt) =2xi{t)-Xi{t-dt)^ai{t)dt^^0{dt^) (A.2) 

is used to update its position. As can be seen, this is a 
third order method, which can be extended also to the 
velocity space |47 



as: 



Vi{t^dt) 



Xi{t) — Xiit — dt) 
dt 



1 



^-dt[llai(t) - 2ai(t - dt)] + O(dt^). 
6 



(A.3) 



The instance when the ith block sticks to the belt is found 
when the relative velocity v^ changes its sign, while the 
instance when the block starts to slip is defined by the 
sign change of Fex — Fgt- A more complicated stochastic 
numerical method was also developed to handle the stick- 
slip dynamics, but was found not to significantly alter the 
presented results. 
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